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Abstract 

We calculate the 7/'-meson energy spectrum in the decay T(IS') — > rj'ggg — > rfX in the 
leading-order perturbative QCD in the static quark limit for the Orthoquarkonium. Our 
principal result is the extraction of parameters of the t]'g*g effective vertex function (EVF) 
involving a virtual and a real gluon from the available data on the hard part of the 7/-meson 
energy spectrum. The perturbative QCD based framework provides a good description 
of the available CLEO data, allowing to constrain the lowest Gegenbauer coefficients B% 
and of the quark-antiquark and gluonic distribution amplitudes of the //-meson. The 
resulting constraints are combined with the existing ones on these coefficients from an 
analysis of the rj — 7 and rf — 7 transition form factors and the requirement of positivity 
of the EVF, yielding S^Oo) = -0.008 ±0.054 and B^Qfi) = 4.6 ±2.5 for = 2 GeV 2 . 
This reduces significantly the current uncertainty on these coefficients. The resulting EFV 
F v 'g* g (p 2 , 0, m^,), including the //-meson mass effects, is presented. 
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1 Introduction 



Experiments involving the production and decays of the 77- and //-mesons are consis- 
tent with the picture that the 77-meson is largely an S77(3)f flavour-octet state, but the 
?/-meson contains a significant amount of a U(3)f flavour-singlet quark-antiquark (qq) 
component in its wave- function [3I3I3II1IS] > in addition, the gluonic (gg) content of the 
?/-meson is substantial This implies that for processes involving gluons and the 7/- 
meson, in particular in the decays of heavy mesons of current phenomenological interest, 
such as b — > sgr] 1 inducing the B — > r]'X s decay |H1 El EMU El EI] and the exclusive 
B -> (r},r}')(K,K*) decays [HII3E3II3II3II3I23I2I] , but also the inclusive decay of the 
Orthoquarkonium T(15) — > gg{g* — > rfg) — > rfX 122], the effective r/-gluon-gluon vertex 
plays an important role. This effective vertex function (EFV) has to be known sufficiently 
well to undertake a quantitative analysis of the data involving the //-meson. Calling this 
EFV F ri > g *g*(pf,pl, m 2 ,), where p\ and p 2 are the virtualities of the two gluons, the QCD 
anomaly provides the normalization of this vertex for on-shell gluons, F v i gg (0, 0, m 2 ,). 
When one or both of the gluons are virtual with relatively large virtuality, the effective 
r}'g*gW vertex can be calculated in perturbative QCD [23I23I23I2EI • 

In this approach, the 7/-meson wave-function is described in terms of two light-cone 
distribution amplitudes (LCDAs) involving the quark-antiquark <jy^){x,Q 2 ) and the glu- 
onic 4>ly\x,Q 2 ) components, where x is the scaled energy of one of the partons of the 
//-meson and Q 2 is the typical hard scale in the vertex. These two components mix 
if the QCD evolution is taken into account. The leading-twist LCDAs of the //-meson 
can be expressed as infinite series of the Gegenbauer polynomials C\ {x — x) (for the 
quark-antiquark) and C^\(x — x) (for the gluonic) components [2T1I2H1I2§1I3 0(13 



y(x,Q 2 ) = 6xx 



1+ MQ 2 )Cl /2 {x-x) 

evenn>2 

$\x,Q 2 )=X 2 X 2 B n{Q 2 )Cn-l{x-x), (2) 

even n>2 

where x — 1 — x, and the following notation is used for the Gegenbauer moments: 

= B^l) (|g|) 7? + & B^l) , (3) 

Bn m = # B^l) ggf + a?>0« (£$f . (4) 

The constrained parameters pn \ pii , 7™ and 7™ are computable and can be found, for 
example, in the Appendix A of Ref. [21]. Usually, one employs an approximate form for 
the //-meson LCDAs in which only the first non-asymptotic terms in both the quark- 
antiquark and gluonic components are kept. Thus, in this approximation 

Q 2 ) = 6^ [1 + 6(1 - hxx) A 2 (Q 2 )} , (jj*\x, Q 2 ) = 5x 2 x 2 (x - x) B 2 (Q 2 ), (5) 

where the explicit forms for C 2 ^ 2 (x — x) and C^ 2 (x — x) have been used. The free param- 
eters B 2 q \/j,l) and £?2 ()Uo), which enter in A 2 (Q 2 ) and B 2 (Q 2 ), are not determined from 
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first principles, and have to be modeled or extracted from a phenomenological analysis of 
data. To that end, a fit of the CLEO and L3 data on the rj — 7 and if — 7 transition form 
factors for Q 2 larger than 2 GeV 2 |33| IHlj was recently undertaken in Ref. [23], yielding 

A 2 (l GeV 2 ) = -0.08 ±0.04, B 2 (l GeV 2 ) = 9 ± 12, (6) 

where the initial scale of the evolution is taken as = 1 GeV 2 . Note, that the coefficients 
A 2 {1 GeV 2 ) and -82(1 GeV 2 ) are strongly correlated. The estimates (jHJ) can be translated 
in terms of the universal free parameters B 2 Q ^ and B 2 9 \ yielding: 

B ( 2 q \l GeV 2 ) = 0.02 ± 0.17, B [ 2 9 \l GeV 2 ) = 9.0 ± 11.5. (7) 

The current determination of these coefficients, in particular B 2 , is rather poor, leading 
to a huge uncertainty in the evaluation of the T]'g*g^ vertex function F ri , gtg (,)(pl,p 2 ,m 2 ,). 

In this paper, we undertake a perturbative-QCD based analysis of the recent data 
on the inclusive process Y(IS') — > rfX in the large 7/-meson energy region published 
recently by the CLEO collaboration [35 , which is expected to be dominated by the 
process Y(IS') — > rj'ggg — > rfX. Moreover, Chen and Kagan [22] have argued that 
the shape of the 7/-meson energy spectrum in this decay is sensitive to the shape of 
the F v > g * g (p 2 ,0,m 2 ,) vertex function, involving the 7/-meson, a real and a virtual gluons 
(see Fig. Q). This sensitivity has already been used by the CLEO collaboration [35] 
to rule out certain models for this vertex, with the CLEO analysis favouring a rapidly 
falling p 2 -dependence of the vertex F v / g * g (p 2 , 0, m 2 ,), in line with the perturbative QCD 
predictions [2*3|l2*4] . Motivated by these observations, we undertake a quantitative analysis 
of the CLEO data to constrain the LCDAs involving the quark- ant iquark and the gluonic 
components of the //-meson. The results of this analysis are presented in terms of the 
coefficients B 2 q \2 GeV 2 ) and B 2 9 \2 GeV 2 ), where we have taken the initial scale as 
^0 = m 5/+Po = ^ GeV 2 , which corresponds to the minimum gluon virtuality p 2 , ~ 1 GeV 2 . 
This analysis is then combined with an earlier analysis of the rj — 7 and 7/ — 7 transition 
form factors j2E] to further constrain the two parameters. As the physical interpretation 
of the function F ri i g * g {p 2 , 0, m 2 ,) is that it represents a probability distribution, much the 
same way as the partonic distributions are in, for example, deep inelastic lepton scattering 
off nucleons, this function must remain positive definite in the entire p 2 range. The 
requirement of positive definiteness of the function F v > g * g (p 2 , 0, m 2 ,) provides additional 

constraints on the parameters B 2 q \2 GeV 2 ) and B 2 9 \2 GeV 2 ), in particular on the latter. 
The combined analysis leads to the following correlated ranges for these coefficients: 

B {q) {2 GeV 2 ) = -0.008 ± 0.054 , B { 2 9) {2 GeV 2 ) = 4.6 ± 2.5 . (8) 

Finally, we use this information to calculate the F v ' g * g (p 2 , 0, m 2 ,) vertex, including the 
Ty'-meson mass effects, relegating the detailed derivation to a subsequent paper jSH] . 

This paper is organized as follows: In section 2, we calculate the branching ratio for 
the process T(15) — > ggg*(g* — > rfg) — > r]'X and the ?/-meson energy spectrum in this 
decay. Numerical analysis of the CLEO data is carried out in section 3, and the resulting 
constraints are combined with the analysis of the ?/— 7 transition form factor to determine 
the LCDAs of the //-meson. Section 4 contains a brief summary. The expressions for the 
matrix element squared for the decay T(15*) — > ggg*(g* — > rfg) are displayed in the 
Appendix. 
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Figure 1: A typical Feynman diagram describing the T(IS') — > ggg*(g* — > t]'g) — > rj'X 
decay. The directed solid line denotes the quark b and antiquark b of the T-meson and 
the double solid line denotes the ?/-meson. 

2 Branching Ratio and ^/-Distribution in the Decay 

T(15) -> rfggg rj'X 

Several processes contribute to the inclusive //-meson production in the T(15')-meson 
decay. The two dominant decays are: T(1S) — > 7* — > qq — > rj'X and T(IS') — > ggg*(g* — > 
rj'g) — > rj'X. The first of these has been estimated from the measured value of the hadronic 
cross-section i?had and the branching ratio B(T(1S) — > yielding B(T(1S) — > 7* — > 

gg) = (8.8 ± 0.3)% [S3], and the CLEO data has been corrected for the qq component. 
In addition, there is also the continuum production from the process e + e~ — >• 7* — > qq — ► 
r^'X, which can be estimated from the e + e~-continuum data below the resonance. This 
data also provides a good profile of the fragmentation in the process 7* — > qq — ► r)'X. 
Typically, the fragmentation processes involving the intermediate qq state populate the 
low- z region, where z = E v r /-Ebeam — 2E 71 r/M is the relative rj'- meson energy expressed in 
terms of the T(lS')-meson mass M. One expects that in the large-z region, the process 
T(1(S) — > ggg*{g* — > grf) dominates. Assuming this, we will concentrate here on the 
intermediate three-gluon {ggg*) state and analyze the -spectrum in the large- z region 
alone. The quality of the fits provides a justification of this procedure a posteriori. 

A typical Feynman diagram describing the decay T(15) — > ggg*(g* — > rfg) — > rj'X is 
presented in Fig. ^ The other 17 diagrams can be obtained from the above one by the 
permutations of the gluons in the intermediate (virtual) and final state. For the matrix 
element calculations, we adopted the static limit for the heavy quark-antiquark pair in 
the orthoquarkonium T(1S), so that both the quark (b) and antiquark (b) carry half of 
the T(15')-meson four-momentum, and the velocity-dependent corrections are neglected. 

The total decay amplitude can be divided into three parts: 

3 

M[T(V) - v'(Pv')9(h)g(k 2 )g(h)} = (9) 

»=i 

where each of the three terms M.i collects the contributions from the diagrams with 
the virtual gluon of the same four- momentum pi = p v > + hi, with p v i and k{ being the 
four-momenta of the ?/-meson and the ith final gluon. The explicit forms of the three 
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contributions Mi follows: 
Mi ■■ 



+ 



+ 
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+ (fi A *f 2 B *)(f3*fv) + (fvf 2 B *)(fi A *ft) + (fvf A *)(ftft) 

' \(fl*ft)(P3fTfvP3) + (fvf 2 B l(P3ftfl*P3) + (fvtf*)M7*ff*P3) 



pi 



+ (fif-nMsfvPs) + (fvfnipsftfiPs) + (fvfnMrfrps) 

where (Iabc (A B,C = 1, . . . , N 2 — 1 are the colours of the gluons) is the symmetrical 
constant of the colour SU(N C ) group with N c = 3; g s is the strong coupling constant; 
(fv)fj,u = 'P^Vu ~ 'PuVv is the polarization tensor of the T(lS')-meson, M, V^, rj^ and ip(0) 
are the mass of the quarkonium, the four-momentum, the polarization vector, and the 
non-relativistic wave-function in the position space evaluated at the origin, respectively; 



(ff 
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ki^i U ki v e i and (/„• 



£ i ivpok p i ef a [i = 1,2,3) are the gluonic field- 



strength tensor and its dual involving the zth gluon with the polarization vector e 
and four- momenta k ifJL \ pi = p v i + ki is the four- momentum of the virtual gluon, and 
p\ = V — pi- In the above equations the notations are as follows: (fif 2 ) = {.fi)pv{.f 2 ) up ', 
(Pif A f 2 Pi) — Pi{fi)pv{f 2 ) up Pi P -, etc. It is worth noting that the matrix element satisfies 
the Bose symmetry under the exchange of gluons, in particular, under the interchange of 
the first and second gluons Mi <-> M 2 while the term Ms remains unchanged. 

The function F ri i g {pf) is the effective vertex function involving the 7/-meson and two 
gluons, one of which is on the mass shell (k 2 = 0). We shall call this interchangeably also 
the if — g transition form factor. A form of this function motivated by the QCD analysis 
of the r]'g*g loop diagram was suggested by Kagan and Petrov j^j: 



F v , g ( P 2 ) = Frf 
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where the function H(p 2 ) was approximated by the constant value H(p 2 ) ~ Hq ~ 
1.7 GeV -1 , extracted from the J/ip — > f]'l decay [7]. In a companion paper |3H| we 
argue that the form ()13J) for the rj' — g transition form factor also emerges in the perturba- 
tive calculations of this function in the hard-scattering approach by keeping the ?/-meson 
mass. In this approach, the dependence of the function H{p 2 ) on the gluon virtuality p 2 
is given by the following expression: 



4^0=) 
m 2 v , 



i + a 2 {q 2 )-\b 2 {q 2 )g { 2 9 \\x) 
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(14) 



where f v > ~ 2f 7S -/\/3 is the //-meson decay constant expressed in terms of the 7r-meson 
decay constant f n ~ 133 MeV, ( = m 2 ,/p 2 , and the function G 2 (1, () has the form 
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ln(l-C). (15) 



In Eq. ((Hj) the scale Q 2 in the strong coupling a s (Q 2 ), and also in the second Gegenbauer 
moments A 2 (Q 2 ) and B 2 (Q 2 ) of the quark-antiquark and gluonic LCDAs of the //-meson, 
is related with the gluon virtuality p 2 , but there is an uncertainty in its precise definition. 
One of the possibilities is to require that the function H{p 2 ) is finite at all values of p 2 
including the singularity point of the EVF (fTHj). p 2 = m 2 ,, which can be done, for example, 
by putting Q 2 = \p 2 \ + m 2 ,. 

The dependence of the function G 2 on the momentum squared p 2 of the virtual gluon 
is presented in Fig. |21 with the value G 2 \l,ti) = 1/6 in the large-|p 2 | asymptotic region. 
Since this result is based on the application of perturbation theory, its validity is restricted 
to the region \p 2 \ > p^, where typically pfy — 1 GeV 2 (or somewhat higher). In view of 
this, we shall set Im G^ = and ignore the structure in Re G^ in the low-|p 2 | region 
seen in Fig. |21 

As already noted, the functions A 2 {Q 2 ) and B 2 {Q 2 ) in Eq. (fH^ are the second Gegen- 
bauer moments of the quark-antiquark and gluonic light-cone distribution amplitudes 
(LCDAs) of the r/-meson. The scale /io appearing in the definitions of these moments 
is set by the initial value in the evolution of the Gegenbauer moments. As the mass of 
the ?/-meson is large, ~ 1 GeV, it is not a good approximation to neglect it. Assum- 
ing further, that the evolution of the Gegenbauer moments is valid for the gluons with 
virtualities larger than 1 GeV 2 , we shall take fxl = Ql ~ 2 GeV 2 . 

The total decay width of the T(15')-meson into the ?/-meson and three gluons, T — > 
ri'ggg, averaged over the quarkonium spin states and summed over the polarizations and 
colours of the final gluons can be written in the form: 

T v , x = T[r^ ri'ggg] = I -L- J- f *1 ^ (16) 



3! (2vr) 8 2M J 2uj x 2u 2 2o; 3 2E V , 



8W(p - kl - k2 - h - Pv ,)lj2\M[r^ ri'gggf 



where the factor 1/3! takes into account the identity of the final gluons. The expression 
for the matrix element squared is rather lengthy and can be found in the Appendix, where 
we have also discussed some technical details of our Monte Carlo integration. 
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Figure 2: The real (solid curve) and imaginary (dashed curve) parts of the func- 
tion GJj (1>C)j where £ = m 2 ,/p 2 , as a function of the momentum squared p 2 of the 
virtual gluon. We use this function only in the \p 2 \ > 1 GeV 2 region. 



The ?/-meson energy distribution function can be defined as follows: 
1 dr n ix{z) 111 1 f dki dk 2 dk 3 dp v > 



dn 
dz 



r <°) dz 



r( 0) 3! (2vr) 8 2M J 2uj x 2uj 2 2cu 3 2E V , 

&9 



(17) 



5 (4) ( p _ h _ h _ h _ _ 2^/M) ~ ^ |A<[T - t/erf 



where is the three-gluon decay width of the T(IS') -meson in the leading order: 



F (0)_ 16 /2 Q \ r o n .3 Mf 2N 1^(0) P 



Here, C F = (iV 2 - 1)/(2JV C ) and S F = (N 2 - 4)/(2JV c ) are the constants of the colour 
SU(N C ) group, and a s (M 2 ) = g 2 (M 2 )/(47r) is the strong coupling estimated at the scale of 
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the T-meson mass. The a s corrections to the decay width Ts g are known since a long time 
and are numerically large [SH]. However, we do not take them into account, as we have 
calculated the decay Y(IS') — > ggg*(g* — > n'g) only in the leading order. One anticipates 
that in the distribution dn/dz, a good part of the explicit 0(az s ) corrections should drop 
out, and we are tacitly assuming that the remaining corrections do not greatly influence 
the energy spectrum derived in the lowest non-trivial order. 

3 Numerical Analysis of the T(15) —> rj'X Data and 
the 7/- Meson LC Distribution Amplitudes 

As the results for the matrix element squared in the decay T(1S) — > ggg*{g* — > Tj'g) are 
not yet available in the literature, we shall give their explicit expressions in this paper. In 
that context we note that it is sufficient to have the expressions for one of the diagonal 
terms \M.i\ 2 (i = 1,2,3) and one of the non-diagonal terms M.iM* (i ^ j), as the 
others can be obtained by using the Bose symmetry. With this, the expressions for the 
components (1/3) \Mi\ 2 , derived from Eq. (fTOjl. and (1/3) Y2MiM%+ c.c, resulting 
from the cross terms in Eq. (fTU|) and (fTTj) in |A4| 2 , are given in the Appendix. 

We start our numerical calculations by reproducing the already known results for the 
//-meson energy spectrum dn/dz |3~H|l2"2"j. as this provides a good consistency check of our 
calculations. For this purpose, the spectrum in Fig. El is calculated with the same set of 
phenomenological parameters as has been used in Ref. [22] for the following three input 
forms 3 : 

a) A slowly falling EVF: F v r g (p 2 ) ~ 2.1 GeV^ 1 [a s (p 2 )/a s (m 2 ,)} where the two-loop 
expression is used for the strong coupling a s . 

b) A rapidly falling EVF of the form p3j1 with the function H(p 2 ) approximated by 
the constant value H{p 2 ) ~ 1.7 GeV -1 . 

c) An intermediate EVF: F v/g (p 2 ) ~ 12.5GeV/(j9 2 + M 2 ) with M = 2.25 GeV. 

The shapes of the 7/-meson energy spectrum resulting from these EVFs are presented in 
Fig. El The experimental points in the decay T(IS') — > n'X shown in this figure are taken 
from Ref. [3E]- Our results are in agreement with the ones in Ref. [22J, and we confirm 
the previous observation that the curve (b) corresponding to the rapidly falling EVF is in 
reasonable agreement with the experimental data in the hard-energy region z > 0.7 [35 . 
Note that the allowed kinematic domain of the variable z is Zq < z < 1 + Zq/4, where 
zq = 2m n ' /M. The region l<z<l + z 2 /4is too small to be resolved experimentally, 
given the statistics of the CLEO data and hence the energy distribution above z — 1 
is included in the last energy bin 0.9 < z < 1.0. 

Let us now turn to the analysis of the perturbative-QCD based form of the effective 
vertex presented in Eqs. (fT3j) - (fT5j) . To fit the parameters (/^o) and B^i/J^) from the 
existing experimental data on the T(1S) — > n'ggg — > n'X decay, it is convenient to first 

3 We thank Alex Kagan for providing us the input parameters used in their analysis of Fig. [3] It should 
be noted that in the rest of our paper we have used the input parameters displayed in Tabled 
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Figure 3: The energy spectrum of the ?/-meson in the decay T(IS') — > rj'ggg — > rj'X 
for the three assumed phenomenological models for the rf g*g vertex given in the text, in 
comparison with the data from the CLEO collaboration [35J. Note that only the data 
points for z = 2£y/M > 0.7 are relevant for comparison with the theoretical models. 



derive an approximate numerical formula for the rj'-meson energy spectrum. This expres- 
sion will come in handy for subsequent comparison with other independent calculations 
and will also facilitate undertaking analysis of the data in future by the experimental 
groups themselves. As the amplitude of the process considered is linear in B!f and B% , 
the energy spectrum dn/dz is quadratic in these parameters. The general form for the en- 
ergy spectrum in the decay T(15) — > rj'ggg — ► rfX can be written as follows (suppressing 
in this equation the scale-dependence of the coefficients for ease of writing): 

fjz, BP, = G.W + C^)B? + C^)B? (19) 

+ C m (z)\Bff + C, 3 (z)BP B« + C sl ,(z)[B!f>} 2 . 

We have generated the theoretical ^/-spectrum, using a high statistics Monte Carlo phase 
space programme, with fixed values of the coefficients B^\ix^) and £?2 (//q), and have 
varied their values over a certain range to numerically determine the dependence of the 
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Table 1: Input parameters used in the numerical analysis. 

M = 9.46 GeV 
m v , = 958 MeV 
/ w = 133 MeV 
ii\ = 2 GeV 2 



Table 2: The coefficients in the interpolating function dh/dz defined in the text for the 
decay Y(IS') — > rj'X. The numbers in the brackets are from the Monte Carlo statistical 
errors. 



z 


C 00 (z) 


C 0q (z) 


Co g (z) 


C m {z) 


C qg (z) 


C 99 (z) 


0.6 -=-0.7 
0.7^0.8 
0.8^0.9 
0.9-=- 1.0 


1.9172(27) 
1.6750(25) 
1.2768(21) 
0.6475(17) 


2.5187(36) 
2.2368(34) 
1.7344(29) 
0.8945(17) 


-0.3108(13) 
-0.2583(12) 
-0.1855(10) 
-0.0866(07) 


0.8744(19) 
0.7953(18) 
0.6343(15) 
0.3396(09) 


-0.1888(07) 
-0.1560(07) 
-0.1097(06) 
-0.0484(05) 


0.0146(06) 
0.0122(06) 
0.0089(05) 
0.0044(05) 



spectrum on these coefficients. Other parameters in our numerical analysis are listed in 
Table d which are the same as the ones used in the analysis of the r)' — 7 and rj — 7 
transition form factors [23], for the sake of consistency, as we are going to make use of 
this analysis. We have not included any errors on m c , m& and Aqq D , as the parametric 
uncertainties from the Gegenbauer coefficients are by far the largest, which we study. The 
coefficients C a b(z) [a,b = 0, q, g] for the four bins of the 7/-meson energy spectrum are 
presented in Table 121 The numbers presented in the parentheses represent the statistical 
error of our Monte Carlo calculations, for which we have used the Monte Carlo phase space 
generator FOWL from the CERN Library of programmes. With the help of the program 
MINUIT [10], the following best fit values of the parameters B^\^q) and B^\/Iq) are 
obtained: 

bPQ® = -0.891{;1, BPOt) = "2.86™, x 2 = 2.45, (20) 
B$»(fi) = -1.09i^, B®{ti) = 11.53-^,. x 2 = 2.37, 

for the last three experimental bins with z > 0.7 (the stated x 2 corresponds to three 
degrees of freedom), and 

B^\^l) = -0.77±8lS, B^\nl) = -4.36^11 X 2 = 24.13, (21) 
B^(nl) = -1.2918:™, B^\^l) = 12.51^:1, x 2 = 23.69, 

for the four experimental data points in the large-z region (z > 0.6). The minimum x 2 of 
the fit in the last case, namely x 2 — 24 for four degrees of freedom, is unacceptably large. 
Thus, as already observed in Ref. [33], only the last three bins in the energy spectrum 
are dominated by the T(IS') — > rfggg —* r]'X decay. Following this, we concentrate on 
the analysis of the last three bins with z > 0.7. The la contours both in the Gegenbauer 
coefficients {B^ (a*o)> -^2 O^o)) anc ^ m ^ ne Gegenbauer moments (A 2 (fJ,l), B 2 (a*q)) are 



m c = 1.3 GeV 
m b = 4.3 GeV 

= 305 MeV 
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Figure 4: The ±1 cr contours (long-dashed curves) in the parameter space B% (/^j) — 
B% (a*o) O e fr frame) and for the Gegenbauer moments B 2 (nl) — A 2 (nl) (right frame) 
resulting from the analysis of the last three experimental bins in the ?/-meson energy 
spectrum in the process T(1S) — > rj'X measured by the CLEO collaboration [SS]- The 
crosses ( x ) represent the two solutions with the minimum x 2 given in the text. The narrow 
solid contours result from the analysis of the data on the rj' — 7 transition form factor, 
scaled from Ref . j2H] . The short-dashed lines in both the frames result by demanding that 
the vertex function F v / g (p 2 ) remains positive for all values of p 2 > m 2 ,, allowing only the 
regions below these lines. Note that all the contours correspond to using fj% = 2 GeV 2 . 

presented as long-dashed curves in Fig. HI The best fit values (|20~j) are denoted by the 
crosses (x). The narrow contours (solid curves) also shown in these figures result from the 
analysis of the data on the ^' — 7 transition form factor [22]. As we are using /Zq = 2 GeV 2 
in our analysis, the results from Ref. [22] have been rescaled from //q = 1 GeV 2 to 
/ip = 2 GeV 2 with the help of Eqs. (j3J) and (|1J) for n — 2. We note that the analysis 
of the 7/ — 7 data, being more sensitive to the quark-antiquark LCDA of the r/'-meson, 
provides a much more precise determination of the parameter B^\/Mq) and hence of the 
Gegenbauer moment A 2 (fil), than what can be determined at present from the CLEO 
data on T(1S) — > rj'X. However, what concerns the coefficient B 2 9 \ij,1), and also the 
Gegenbauer moment B 2 (n1), the CLEO data on T(1S) v'X, despite its statistical 
limitations, has cut out some of the allowed parameter region from the rj' — 7 analysis. 
The additional constraint that follows from demanding that the vertex function F v / g (p 2 ) 
remains positive definite for all values of p 2 > m 2 , is shown through the short-dashed lines 
in these figures, admitting only the parameter space below these lines. We will discuss 
this constraint in more detail below. 

A blow up of the overlapping region in the parameter space (5^ (/Zq),^ (a*o)) ^ s 
shown in Fig. 03 Here, the la contours following from the rj' — 7 transition form factor 
analysis and the T(IS') — > rj'X decay fit are shown through the solid and dashed curves, 
respectively. In drawing the allowed parameter space from the T(1S) — > rj'X decay fit, we 
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Figure 5: The combined fit for the parameters B 2 (a*o = 2 GeV 2 ) and B% (/i 2 , — 2 GeV 2 ) 
from the data on the rj — 7 transition form factor (solid curve) and T(IS') — > r/X decay 
(long-dashed and short-dashed curves) with the requirement that the vertex F v > g (p 2 ) re- 
mains positive definite in the entire p 2 region. The resulting la contour (combined best 
fit) is shown by the yellow (shaded) region. The point denoted by (^) corresponds to the 
asymptotic light-cone distribution amplitude. 



have imposed the additional condition that the r)'g*g vertex function F v > g (p 2 ) for p 2 > m 2 , 
has the same sign as the corresponding function calculated with the asymptotic forms of 
the r/-meson quark-antiquark and gluonic LCDAs which is defined as positive-definite in 
the entire p 2 region. With Eq. (|13|) taken into account, this condition implies H(p 2 ) > 
and results [with the help of the explicit from (fTfj) of the function H(p 2 )] in the following 
inequalities: 

S6 

B 2 (A*o) < Y [l + A 2 ^D], (22) 

{g) , 36 + (36 -5^)5^^) 
B 2 (ji ) < 



5 - 36p { 2 9) 

where the function G 2 9 \l, C) ^Hfy is approximated by its value 1/6 in the large- \p 2 \ asymp- 
totics. The values p [ 2 q) = 2.86 and p { 2 9) = -0.01 were taken for the constrained parameters 
in the numerical analysis. The positivity constraint removes the larger values of B^ip 2 ,) 
above the short-dashed curve in Fig. which would otherwise force the vertex function 
to cross the zero at some value of p 2 and become negative. This is exemplified below 
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Figure 6: The ?/-meson quark- ant iquark and gluonic LCDAs 4>3 ( function 
of x (left frame), and the resulting rf — g transition form factor F v / g (p 2 ) (right frame). 
The solid curves correspond to the asymptotic quark-antiquark LCDA, while the LCDAs 
for the combined best fit values of the Gegenbauer coefficients given in Eq. ()23|) are drawn 
as dotted curves. The LCDAs with the values B% = 0.15 and = 13.5, which are 
allowed within la from the analysis of the data on the rf — 7 transition form factor , 
are presented as the dashed curves. Note, that for this case, the function F v / g (p 2 ) is no 
longer positive-definite, as shown in the right frame. 

for specific values of the Gegenbauer coefficients. The resulting combined best fit of the 
parameters -63 (//q) and B% (a*o) ^ s snown as the coloured (shaded) region. This yields 
the following correlated values: 

B ( 2 q) = -0.008 ± 0.054, B { ° ] (fx 2 ) = 4.6 ± 2.5, (23) 
A 2 {f4) = -0.054 ± 0.029, # 2 (//<!) = 4 - 6 ± 2 - 7 > 

with the central values having x 2 — 2.66 for three degrees of freedom. For comparison, 
the point shown as (^) , corresponding to the asymptotic LCDA, has \ 2 = 8.41 for three 
degrees of freedom. Thus, the asymptotic form of the r/-meson LCDA is not favoured by 
our analysis. 

The shapes of the quark-antiquark and gluonic LCDAs are presented in Fig. |H1 (left 
frame); the resulting rf — g transition form factor F v / g (p 2 ) corresponding to these LCDAs 
is also shown in this figure (right frame). The solid and dotted curves in this figure 
correspond to the asymptotic LCDA and the combined best fit values of the Gegenbauer 
coefficients (|23jh respectively. We also show, for the sake of illustration, representative 
LCDAs with the Gegenbauer coefficients having the values B^(2 GeV 2 ) = 0.15 and 
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Figure 7: The rf—g transition form factor F v / g (p 2 ), defined in Eq. (fT3*j) (left frame), and the 
function m 2 ,H(p 2 ), defined in Eq. (jl4j) (right frame), as functions of the gluon virtuality p 2 . 

The various curves correspond to the input values of the Gegenbauer coefficients B^ (/^q) 
and B^\jjL^) resulting from the combined best fit shown in Fig. The solid curves 
correspond to the central values of the fit, while the dashed curves are drawn for the 
parameters corresponding to the uppermost and lowermost allowed values of the combined 
best fit contour shown in Fig. |SJ 

B^\2 GeV 2 ) = 13.5. These parametric values are taken from the analysis in Ref. |25] . 
but are in conflict with the positive definiteness of the vertex function F n i g (p 2 ), as shown 
by the dashed curve for F v i g (p 2 ). In fact, the constraint of positivity on the effective vertex 
function cuts out regions in the parameter space b\ 9 '(2 GeV 2 ) > 0.045 and B^^ GeV 2 ) > 
7.1, as otherwise the contribution from the gluonic LCDA starts to dominate, which 
makes the vertex function F v / g (p 2 ) negative for some range of p 2 considered here and in 
Ref. [25J. Thus, positivity criterion provides an effective constraint on the magnitude of 
the coefficient B^ (/J^), reducing significantly the otherwise allowed range in Eq. ((7|). In 
the context of our analysis, we note that the asymptotic //-meson LCDA provides a fair 
description (though not the best fit) of the current data on T(15) — > ggg*(g* — > rj'g) — > 
g'X, and hence one anticipates that the subasymptotic corrections in the LCDAs, and 
the vertex function F n i g (p 2 ), while important in the analysis of data, are not dominant. 4 
The EVF F v ' g (p 2 ) and the function m 2 ,H(p 2 ) for the combined best fit values (J23*j) 
are presented in Fig. The //-meson energy spectrum for the EVF motivated by the 
perturbative QCD analysis is presented in Fig. |H1 in comparison with the experimental 
data [35] and the spectrum corresponding to the "rapidly falling" form of the EVF [22] 
labeled by the constant value of the function H(p 2 ) = H ~ 1.7 GeV -1 . The second dotted 
line labeled by if as (p 2 ) corresponds to the EVF when only the asymptotic form of the r]'- 

The approximate validity of the asymptotic transition form factors involving the ir-, rj- and 7/-mesons 
compared to the data, observed in Refs. |37I38| . also suggests that the coefficient ((J-q) in the //-meson 
is bounded by this data. The gluonic components of the ir- and the ?7-meson LCDAs are small in any 
case. 
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meson LCDA {B^ = and B<f = 0) is taken into account. The yellow (shaded) region 
demarcates the spectrum with the Gegenbauer coefficients having values in the range of 
the combined best fit (j2SJ)- The blue (solid) curve lying inside this region corresponds to 
the best fit values (J2Uj) obtained from the analysis of the Y(IS') — > rj'X decay only. From 
Fig. |S1 it is seen that the rapidly falling phenomenological EVF with H = 1.7 GeV -1 
gives a harder //-meson energy spectrum for the large energy z region compared to the 
CLEO data. The spectrum with the asymptotic form of the ?/-meson LCDA is well 
correlated with the experimental point in the bin 0.7 < z < 0.8, but overestimates the 
data in the other two bins. It should be noted that the last two bins (especially the 
bin 0.8 < z < 0.9) are very uncertain in the current data, which has to be statistically 
improved to draw sharper conclusions. 

0.015 



o.oio- 
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z 



Figure 8: Energy spectrum of the r/'-meson in the decay T(1S) —>■ rj'X. The upper dotted 
curve corresponds to a constant value of the function H(p 2 ) — Hq ~ 1.7 GeV -1 , and the 
curve marked as if as (p 2 ) corresponds to the asymptotic form of the ?/-meson LCDA (i.e., 
£?2 = and = 0). The spectrum with the Gegenbauer coefficients in the combined 
best-fit range of these parameters is shown by the yellow (shaded) region. The blue (solid) 
curve inside this region corresponds to the best fit values of the Gegenbauer coefficients 
from the analysis of the T(IS') — > rj'X data alone. 
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4 Summary 



We have calculated the r/'-meson energy spectrum in the decay T(1S) — > ifggg — > tj'X 
in leading order perturbative QCD in the static quark limit for the Orthoquarkonium. 
Assuming some phenomenological vertex functions, our results are in agreement with the 
ones obtained earlier in Ref. [221- I n our analysis, the ?/-meson is described in the leading- 
twist (twist-two) accuracy in terms of the quark-antiquark and gluonic LCDAs for which 
the asymptotic and the first non-asymptotic terms are taken into account. In this ap- 
proximation, the T]'g*g EVF depends essentially on the Gegenbauer coefficients B^\hq) 
and B^\nl). They are determined from the CLEO data on T(1<S) — > rfX in the large-;? 
region (z > 0.7) of the ?/-Hieson energy spectrum [HE], which is well explained by our 
perturbative QCD approach. However, the resulting la contour in the Gegenbauer coef- 
ficients have a large dispersion. Combining the analysis of the T(1S) — > rj'X data with 
an earlier analysis of the (7/, 77) — 7 transition from factors in a similar theoretical frame- 
work j2Sj, and requiring additionally that the vertex function F v ' g (p 2 ) remains positive 
definite in the entire p 2 > m 2 , region, yield much improved determination of the Gegen- 
bauer coefficients, yielding B { 2 q) (2 GeV 2 ) = -0.008 ± 0.054 and B { 2 9) (2 GeV 2 ) = 4.6 ± 2.5. 
Our analysis improves the phenomenological profile of the LCDAs in the ry'-meson, and in 
turn yields a better determination of the vertex function F v / g (p 2 ) compared to the earlier 
estimates of the same. The resulting function (the rf - gluon transition form factor) is 
presented in this paper including the r/'-meson mass effects. 
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A Matrix Element Squared for the Decay T(1S) — > 
9 99* (9* rfg) 

In Sec. El we have displayed explicitly the terms Aii (i = 1,2,3) in the decomposition of 
the total amplitude of the decay T(15) — > r]'ggg in the leading-order perturbative QCD 
[Eqs. (TTU Jl - ffT^j) ]. Squaring this amplitude results in very lengthy expressions. Neverthe- 
less, it is possible to present compact analytical expressions for some parts of the matrix 
element squared. In particular, the square of the term Aii, summed over the polarizations 
and colours of the final gluons and averaged over the polarization states of the T-meson, 
can be written as follows: 



— CpBpTT a s (M ) \F n , g { Px )\ (pA;2)2( p fc3)2to)2(p?)2 ' (24) 



15 



where the well-known property of the symmetrical constants (Iabc of the colour SU (N c )- 
group has been used: 



£ AT-i A 

^ d-ACD dj3CD = ~tt — 5ab = 2B f 5ab, (25) 



C,D=1 ' c 



and the dynamical function Ju(V, k^Pi) has the form: 

Jn(V,h, Pi ) = AM^hp^^ihhYiVp.y + ihp^iVhr + ihp^iVh) 2 ] (26) 

+ AM 2 {p\) 2 {k 2 k 3 ) 2 [(hk 2 ) 2 + (hk 3 ) 2 ] 

+ 8M 2 pj(p 1 k 1 ) 2 (k 2 k 3 ) [(k 2 k 3 ) 2 - (k 2Pl )(h Pl )] 

- SMYMh^h) 2 [(hk^hp,) + (hk^ihp,)] 

+ 2M 2 pj(M 2 -pj) {[(hk 2 )(k 3Pl ) - (hk 3 )(k 2Pl )] 2 - ( Pl h) 2 (k 2 k 3 ) 2 } 

+ 2p\{M 2 -pl) [M 2 {k 2 k 3 ] + 2{Vk 2 ){Vk 3 )] 

x [pKhk^ik^) - (fcipi)(fciA;2)(A;aPi) - (k 1 p 1 )(k 1 k 3 )(k 2 pi)] . 

It is easy to check that this function is invariant under the interchange (k 2 ,p 2 ) <-> (& 3 ,p 3 ), 
which reflects the Bose symmetry for the set of the Feynman diagrams considered. The 
other two dynamical functions J 22 (V, ki,pi) and J 33 (V, ki,pi) originating from the squares 
of the terms Ai 2 and Ai 3 , respectively, can be obtained from Jn(V, ki,pi) by the obvious 
replacements: (ki,pi) <-> (k 2 ,p 2 ) and (ki,pi) <-> (& 3 ,p 3 ). 

The differential width of the 1 — > 4 decay has a non-trivial dependence on five variables 
(three angles can be integrated out trivially) in the rest frame of the decaying particle. 
In the Monte Carlo generator, the matrix element squared is rewritten in terms of the 
following dimensionless variables: 



_ 2{Vh) p 2 

which satisfy the relation: 



*« = ^. * = W> * = 1 ' 2 ' 3 ' ( 2 7) 



x 1 + x 2 + x 3 + y x + y 2 + y 3 = 2 + z 2 /4, (28) 

where z = 2m rj i jM is the relative mass of the r/'-meson. 

The experimentally measured quantity is the r/'-meson energy spectrum - the number 
of T — > rj'X events Arii inside the energy bin with the central value Zi and the width Az. 
This can be calculated theoretically in the Monte Carlo approach, which we have used. 
In particular, as an input for the calculation of the contribution coming from |.Mi| 2 , the 
following equation was used: 



Az 

Zi-Az/2 



\MF^ g ( yi )\ 2 Fujxuyufj.) 

6^(n 2 - 9)M4 [ yi x 2 x 3 (2 - 2 Vl - x 2 - x 3 )} 2 ' 1 ) 
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where the representation (jUj) for |*Mi| 2 was put into the definition of dn/dz (|T7|). The 
dimensionless functions T\\{x^y^\i) is related to the functions Ju(V, ki,pi) as follows: 

M 14 

Ju{V, ki,pi) = -j^Fii[Xi,yi,lj)- (30) 

The first integral in Eq. (|29|) serves for a selection of the generated Monte Carlo events 
with the 7/-meson energy being inside the bin Z{ — Az/2 < z < Z{ + Az/2. 

In the same way all the other contributions from the matrix element squared were 
implemented in the numerical analysis. In particular, the contribution to the energy 
spectrum coming from the product (1/3) M\M. 2 + c.c. can be written as: 



Anj 12) 
Az 



Zi+Az/2 



[ dz6(z-2E v ,/M) f^i^^l^Ls^iP-h-h-ks-p.) 
J K v ' > J 2cu! 2co 2 2cu 3 2E V , y 123^; 

M 2 F v , g (y 1 )F r ,, g (y 2 ) 2J zr i 2 (x^, /i) 

647T 5 (7r 2 - 9)M 4 y 1 y 2 x 1 x 2 x 2 3 (2-2y 1 -x 2 -x 3 )(2-2y 2 -x l -x 3 y 1 ' 

However, we have not been able to find compact expressions for the non-diagonal terms 
(1/3) ^2 A4iA4* (i 7^ j). We present here instead the dimensionless function T\ 2 {xi^y^\x) 
in a form of a series in powers of z^/4 = m 2 ,/M 2 : 

v2\ k 



/ 2 \ K 

Fi 2 (x l ,y l ,z )=J2fi2(^y l ) [^) , (32) 

fc=0 ^ ' 

where functions /f 2 ( follows: 



A° 2 = i {(Ax 12 ) 4 (A z/12 ) 2 -4(Ax 12 ) 3 (A z/12 ) 3 + 4(Ax 12 ) 3 A z/12Z/3 2 (33) 

+ 8(Ax 12 ) 3 Ay 12 y 3 (x 12 + x 3 - 2) 

+ 4(Ax 12 ) 3 Ay 12 (x 2 2 + 2x 12 x 3 - 4x 12 + x 2 - 4x 3 + 4) 

+ 5(Ax 12 ) 2 (Ay 12 ) 4 + 2(Ax 12 ) 2 (Ay 12 ) 2 y 3 (-2x 12 - 4x 3 + 7) 

+ 2(Axi 2 ) 2 (Ay 12 ) 2 (-x 2 2 - 6x12X3 + 7xi 2 - 2x 2 + 9x 3 - 8) 

- 2Ax 12 (Ay 12 ) 5 + 4Ax 12 (Ay 12 ) 3 y 2 (-x 3 - 2) 

+ 4Axi 2 (Ayi 2 ) 3 ?/ 3 (-xi 2 x 3 - 2xi 2 + 2x 3 ) + 4Axi 2 (Ay 12 ) 3 xi 2 (-xi 2 + x 3 + 1) 

+ 2Axi 2 Ay 12 ?/ 4 (2x 3 + 5) + 4Ax 12 Ay 12 ?/ 3 (3xi 2 x 3 + 2x i2 + 2x 2 - 12x 3 + 12) 

+ 4Ax i2 A|/i 2 ?/|(3x 2 2 x 3 + 4x i2 x 2 - 17x i2 x 3 + 19x i2 + x 3 - llx 2 + 36x 3 - 36) 

+ 4Axi 2 Ayi 2 ?/ 3 (x? 2 x 3 + 2x 2 2 x 3 - 12x 2 2 x 3 + 14x 2 2 + x X2 x\ - 16xi 2 x 3 

+54x i2 x 3 - 52x i2 - 4x 3 + 28x 3 - 64x 3 + 48) 

+ 2Axi 2 Ay 12 (— x 4 2 — 6x 3 2 x 3 + 14x 3 2 — 8x 2 2 x 2 + 44x 2 2 x 3 — 56x 2 2 — 2x\ 2 x\ 

+30xi 2 x^ - 96xi 2 x 3 + 88x12 + x\ - 24x^ + 64x 3 - 48) 

+ 4(At/ 12 ) 4 2 /3 2 + 2(A 2 / 12 ) 4 2/3 (3x 12 -l) 

+ (Ay 12 ) 4 (3x 2 2 + 2x i2 x 3 - 6x i2 + 2x 2 - 6x 3 + 4) 
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- 8(Ay 12 ) 2 y* + 4(Ay 12 ) 2 ^( a;i2a ;3 _ 4xi 2 + 2x 2 + 2x 3 - 7) 

+ 4(A?/i2) 2 ?/f (2x 2 2 x 3 — 2x 2 2 + 3xi 2 x 2 — 4xi 2 x 3 — 7x 12 — 9x 2 — 3x 3 + 20) 

+ 2(At/i 2 ) 2 2/ 3 (2x 3 2 x 3 - 2x\ 2 + 2x\ 2 x\ - 10x 2 2 x 3 - 7x 2 2 - 20xi 2 x 2 - 4xi 2 x 3 

+48xi2 + 10x 2 + 32x 3 - 52) 
+ (Ayi 2 ) 2 (x\ 2 — 14x 3 2 — 12x 2 2 x 3 — 10x 2 2 x 3 + 56x 2 2 — 8xi 2 x 3 + 36xi 2 x 3 

+32xi 2 x 3 - 88x12 - 4x^ + 2Qx\ - 32x 2 - 24x 3 + 48) 

(Ax 12 )Y 3 + 2(Ax 12 ) 4 y 3 (-x 12 - x 3 + 2) 
+ (Axi 2 ) 4 (-x 2 2 - 2xi 2 x 3 + 4xi 2 - x\ + 4x 3 - 4) 
+ 7(Axi 2 ) 2 2/ 3 4 + 2(Axi 2 ) 2 y 3 3 (8xi2 + 10x 3 - 19) 
+ 2(Axi 2 ) 2 ?/I(6x 2 2 + 20xi 2 x 3 - 33xi2 + llx 2 - 43x 3 + 40) 
+ 2(Axi 2 ) 2 t/ 3 (2x? 2 + 16x 2 2 x 3 - 21x 2 2 + 20xi 2 x^ - 70xi 2 x 3 + 56xi 2 + 6x 3 

-37x 2 + 72x 3 - 44) 
+ Ax 2 2 (x 4 2 + 12x^ 2 x 3 - 14x 3 2 + 24x 2 2 x 2 - 78x 2 2 x 3 + 56x 2 2 + lQx l2 x 3 3 

-90xi 2 x 2 + 160xi 2 x 3 - 88x12 + 3x 3 - 26x 3 + 80x 3 - 104x 3 + 48) 
+ Ayl + 2^(-2xi 2 x 3 + 5xi2 - 4x 2 - 4x 3 + 15) 

+ yj(— 16x 2 2 x 3 + 9x 2 2 — 36xi 2 x 2 + 38xi 2 x 3 + 50xi 2 — I6X3 + 82xg — 14x 3 — 
+ 2yf (-12x 3 2 x 3 + 2x\ 2 - 30x 2 2 x 3 + 56x 2 2 x 3 + 15x 2 2 - 22x X2 x\ + 126xi 2 x 2 

-68xi 2 x 3 - 64xi2 - 4xg + 52x^ - 126x 2 + 24x 3 + 52) 
+ y 2 (-16x 4 2 x 3 + x\ 2 - 44x? 2 x 3 + 104xf 2 x 3 + lOx 3 ^ - 40x 2 2 x 3 ; + 298x 2 2 x 2 . 

-250x 2 2 x 3 - 64x 2 2 - 12x12X3 + 228x12X3* - 644xi 2 x 2 + 288xi 2 x 3 + 104xi 2 

+40x 4 - 260x 3 + 456x 2 - 152x 3 - 48) 
+ 2y 3 x 3 (— 2x\ 2 — 6x 4 2 x 3 + 19x 4 2 — 6x 3 2 x 2 + 66x 3 2 x 3 — 74x 3 2 — 2x 2 2 x 3 

+68x 2 2 x 2 - 227x 2 2 x 3 + 144x 2 2 + 21xi 2 x 3 - 156xi 2 x 2 + 304xi 2 x 3 

-136xi2 - 17x 3 3 + 88x 2 - 132x 3 + 48) 
+ a: 3 (7x 4 2 + 16x 3 2 x 3 — 46x 3 2 + 13x 2 2 x 2 — 78x 2 2 x 3 + 112x 2 2 + 6xi 2 x 3 

-46x i2 x 2 + 128xi 2 x 3 - 120xi2 + 2x 4 - 14x 3 + 44x 2 - 72x 3 + 48)}, 

= - j - 8(Axi 2 ) 3 Ai/i2i/ 3 + 4(Axi 2 ) 3 Ai/i2(-xi2 - x 3 + 2) 

+ 2(Axi2) 2 (Ayi 2 ) 2 y 3 (x 3 + 7) + (Axi 2 ) 2 (Ayi 2 ) 2 (2xi2X 3 + 4x i2 - 2x 3 - 9) 

- 6Axi 2 (Ayi 2 ) 3 y 3 + 2Axi 2 (Ayi 2 ) 3 (x 3 + 2) - 2Axi 2 Ayi 2 ^(2x 3 + 3) 
+ 2Axi 2 Ai/i2i/ 3 (-4xi2X 3 - 2xi2 - 2x 2 + 7x 3 - 10) 

+ 2Axi 2 A|/i2i/ 3 (-2x 2 2 x 3 - 3x 2 2 - 2xi 2 x 3 + 2xi 2 x 3 + 2x i2 - x 2 + 8x 3 ) 
+ 2Axi2 Aj/12 (-2x 3 2 + ^12^3 + 8x 2 2 + 4xi 2 x 2 - 2x i2 x 3 - 12x i2 + x 3 
-2x 2 - 4x 3 + 8) 

- (At/i 2 ) 4 (xi2 + 1) + 8(Ay 12 ) 2 yl + 2(At/i2) 2 y 3 2 (-xi 2 x 3 + 9x i2 + 2x 3 + 5) 
+ 4(Ayi2) 2 y 3 (-x 2 2 x 3 + 2x 2 2 + 2x i2 x 3 - 2xi2 - 6x 3 ) 

+ (Ayi 2 ) 2 (-2x 3 2 x 3 + 4x 3 2 + 6x 2 2 x 3 - 15x 2 2 - 6xi 2 x 2 + 16xi 2 - 4x 3 
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+18x 2 3 - 8x 3 - 4) 
+ (Ax 12 ) 4 y 3 + (Ax 12 ) 4 (x 12 + x 3 - 2) + 2(Ax 12 ) 2 y 3 3 (-x 3 - 8) 
+ (Axi 2 ) 2 ^(-6a; 12 a;3 - 22x i2 - Ax 2 3 + 4x 3 + 27) 

+ 2(Axi 2 ) 2 y 3 (— 3xf 2 x 3 — 3x 2 2 — 4xi2X 2 + Ax\ 2 x 3 + \\xi 2 — x\ + 6x3 — 3x 3 — 8) 
+ (Axi 2 ) 2 (— 2xf 2 ^3 — 4x 2 2 x 2 + 10x 2 2 x 3 — x\ 2 — 2x12X3 + 18x12X3 — 30xi 2 x 3 

+8x12 + 8x\ - 33x 2 + 40x 3 - 12) 
- 8y\ + yl(6x 12 x 3 - 21x 12 + 8x 2 + 4x 3 - 41) 

+ 2yf (10x 2 2 x 3 - 9x 2 2 + 14xi 2 x 2 - 24xi 2 x 3 - 24x i2 + 4xg - 24x 2 + 12x 3 + 44) 
+ 2/3(24x13X3 - 6xf 2 + 36x 2 2 x 2 - 88x 2 2 x 3 - 19x 2 2 + 14x12X3 - 94xi 2 x 2 

+56x12X3 + 112xi2 - 16x3 s + 50x i + 32x s ~ 1°°) 
+ y 3 (12x A 2 x 3 — x\ 2 + 20x^ 2 x 2 — 56x^x3 — 10xf 2 + 8x 2 2 Xg — 86x 2 2 x 2 + 82x 2 2 X3 

+64x 2 2 - 32x12X3 + 152xi 2 x 2 - 56x12X3 - 104xi2 + 40xg - 120x 2 + 40x 3 + 48) 
+ x 3 (2xi 2 + 4x12^3 - 19x^2 + 2xi 2 a; 2 - 34xi 2 a;3 + 74x^ 2 - 8x 2 2 x 2 + 97x 2 2 x 3 - 144x 2 



+\\x l2 x\ - 8x l2 x\ - 104xi 2 x 3 + 136xi2 + 4x^ - 25x^ + 32x^ + 28x 3 - 48) 



+ 4Axi2At/i2t/|(x 3 - 2) + 4Axi 2 Ai/i22/3 (^12^3 - 2x i2 - 6x 3 + 8) 

+ 4Ax 12 Ay 1 2(3x 2 2 - x 12 x 3 - X12 + 4x 3 - 4) - 8{Ay X2 ) 2 y\ - 12(A2/i 2 ) 2 t/3(2a;i2 + 1) 

+ 4(A2/i 2 ) 2 (-2x 2 2 + x 12 - x 2 + 5x 3 ) + 3(Ax 12 ) 4 + 2(Ax 12 ) 2 y3 2 (4x 3 + 23) 

+ 2(Ax 12 ) 2 y 3 (8a;i2X3 + 20xi 2 + 4x 2 - 6x 3 - 17) 

+ 2(Axi 2 ) 2 (4x 2 2 x 3 + 4X12X 2 ; + 4xi 2 x 3 - 9x i2 + 5x 2 - 27x 3 + 16) 

+ 2Ay\ + 4yf (-3x12X3 + 17xi2 - 2x 2 + 2x 3 + 21) 

+ 22/|(-16x 2 2 x 3 + 27x 2 2 - 10xi 2 x 2 + 46xi 2 x 3 + 22x 12 + 22x 2 - 50x 3 - 52) 

+ 2y 3 (— 14x 12 X3 + 6x 12 — 10x 12 x 3 + 44x 12 x 3 — 7x 12 + I8X12X3 — 12x12X3 

-32xi2 + 10x 2 3 - 56x 3 + 44) 

— 8xi 2 x 3 + x\ 2 — 8xi 2 x 2 + 8x12^3 + 10x^2 + 18x 2 2 x 2 + 86x 2 2 x 3 — 64x 2 2 

+ 28X12X3 1 - 52x 12 x 2 - 176xi 2 x 3 + 104xi2 + 12x 4 - 68x 3 3 + 88x 2 + 72x 3 - 48), 



- 2(Axi 2 ) 2 2/ 3 (a:3 + 6) + (Ax 12 ) 2 (-2x 12 x 3 - 6x 12 + 8x 3 + 1) 

- 8yf + 2y 2 3 (x 12 x 3 - 13xi 2 - 2x 3 - 11) + 2y 3 (2x 2 2 X3 - 9x 2 2 - 8x 12 x 3 + 20x 3 + 4) 

+ 2xf 2 x 3 - 2x\ 2 - \2x\ 2 x 3 + llx 2 2 + 6x12X3 + 24xi 2 x 3 - 16xi 2 + 4xg - 22x 2 + 4 




(35) 




(36) 




(37) 
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fu = -i(*i3 + l), (38) 

where X12 — x\ + £2, Axi2 = 21 — 22, and A?/i 2 = y\ — yi. It is easy to see that the 
equations for f* 2 are symmetric under the interchange (x\, y±) (x%, y<i). 

The expressions for the other Tij{x^ yi, fi) can be derived from the ones given above 
by using the Bose symmetry. 
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